Universal scaling at field-induced magnetic phase transitions 
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We study field-induced magnetic order in cubic lattices of dimers with antiferromagnetic Heisen- 
berg interactions. The thermal critical exponents at the quantum phase transition from a spin liquid 
to a magnetically ordered phase are determined from Stochastic Series Expansion Quantum Monte 
Carlo simulations. These exponents are independent of the interdimer coupling ratios, and converge 
to the value obtained by considering the transition as a Bose-Einstein condensation of magnons, 
Qbec = 1-5. The scaling results are of direct relevance to the spin-dimer systems TlCuCls and 
KCuCla, and explain the broad range of exponents reported for field- induced ordering transitions. 
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A quantum spin liquid is characterized by a finite mag- 
netic correlation length, which is inversely proportional 
to the energy gap, Aq, between the singlet ground state 
and the lowest triplet excitation. Application of an ex- 
ternal magnetic field causes a linear reduction of Aq by 
the Zeeman effect. At zero temperature there is a field- 
induced quantum phase transition (QPT), at a lower 
critical field hci = ~ Aq, from the spin liq- 

uid to a magnetically ordered state. In this phase the 
uniform magnetic polarization then rises continuously 
from zero to saturation at a second critical field hc2- 
Such transitions have been observed in a variety of quan- 
tum antiferromagnets, including Cu2(C5Hi2N2)2Cl40, 
KCuClall 0, and TlCuClsH 11 0, which consist of 
weakly coupled, low-dimensional structural units such as 
spin ladders or dimers. In the ordered phase, the uniform 
magnetization is accompanied by a long-ranged antiferro- 
magnetic (AF) spin ordering in the plane perpendicular 
to the field, which persists up to a finite transition tem- 
perature, Tc, determined by the smallest coupling con- 
stant between the magnetic structural units. 

The compounds KCuCla and TlCuCls provide an ex- 
cellent illustration of the field-induced QPT from a dimer 
spin liquid to a three-dimensional (3D) ordered phase. 
These systems consist of dimer units formed by pairs 
of S" — 1/2 Cu^"*" ions with significant interdimer su- 
;erexchange interactions in all three spatial directions 
0| . Because these interactions are considerably weaker 
in KCuCla than in TlCuCls, the magnon modes of the 
former have a narrow band and a large spin gap (Aq = 
2.6meV), while the latter represents a wide-band system 
with a small gap (Aq = 0.7meV). The wave vector Q 
of the 3D AF order is the location in momentum space 
of the band minimum, Aq, of the triplet magnon exci- 
tations, which are gapped in the disordered phase. The 
field-induced QPT may be regarded as a condensation 
of magnons in the lowest triplet branch, which becomes 
the gapless excitation, or Goldstone mode, correspond- 
ing to rotations of the staggered magnetic moment about 
the field axis and reflecting the broken symmetry of the 



system jg. The continuous evolution of the ordered mo- 
ments is then equivalent to a field-dependent density of 
condensed magnons, and the magnetic field represents an 
experimentally accessible tuning parameter for the QPT. 

The suggestion that this transition may be described 
as the Bose-Einstein condensation (BEG) of hard-core 
bosons 7] has motivated extensive experimental and the- 
oretical studies of the critical properties of such systems. 
Particular attention is given to the scaling exponent a 
which determines the power-law dependence 
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of the ordering temperature on the field h in the vicinity 
of the QPTs at he = hd and hc2- Different experimental 
studies indicate that the precise value of this exponent at 
hci depends on the material, as well as on the measuring 
technique. The exponents a « 1.8 — 2.3 extracted from 
experimental studies of KCuCla and TlCuClg '5,11111 
are all larger than the value aeEC = 3/2 obtained from 
the BEG scenario for 3D systems in which the magnon 
modes have quadratic dispersion at the band minimum 
FP]. The analytic bond-operator approach, applied with 
the parameters of KCuCls and TlCuCla, also yields the 
respective values 2.6 and 1.8 0, while recent Quantum 
Monte Carlo simulations for weakly coupled dimers 
suggest that a ^ 2.3. In both experiments and simula- 
tions, precise extraction of the critical exponent is com- 
plicated by the fact that typically only rather few data 
points are available close to the QPT. In addition, the ex- 
ponent may not refiect universal critical behavior if the 
temperature exceeds an (unknown) upper bound for the 
validity of the BEC approach. However, it is also possi- 
ble that the scaling behavior is simply not universal, in 
that the exponent depends on certain system parameters, 
most notably the strength of the interdimer interactions. 

To analyze the scaling behavior associated with field- 
induced magnetic order, we study a Heisenberg model of 
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FIG. 1: Phase diagram of 3D AF coupled dimers in a mag- 
netic field. The sohd line represents /ici and the dot-dashed 
line /ic2. Circles and triangles are zero-temperature fits ob- 
tained from SSE QMC results. Inset: one layer of coupled 
dimers. 
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FIG. 2: Scaled staggered magnetization at constant field 
h = 0.8J and interdimer coupling J' = 0.25J, for a range 
of system sizes. Tc is determined from the intersection of 
the four curves (main panel), and its error from the extremal 
intersection points (inset). 



coupled spin- 1/2 dimers with the Hamiltonian 

// = ^J.,S,.S,-/i^5f. (2) 

The AF coupling constants J and J' between nearest- 
neighbor sites i and j denote respectively the intra- and 
interdimer superexchange interactions for the cubic ge- 
ometry of which a planar section is illustrated in the inset 
of Fig. n Because the K /TlCuCls structure has no mag- 
netic frustration, the AF cubic dimer system is expected 
to illustrate the same generic similarities and differences 
between large- and small-gap systems as those obtained 
in the experimental studies to date. 

We take as the order parameter the component of stag- 
gered magnetization perpendicular to the applied field, 

= (3) 

where L denotes the linear system size and the per- 
pendicular staggered structure factor, defined as 

This quantity is calculated using the Stochastic Series 
Expansion Quantum Monte Carlo (SSE QMC) method 
of Refs. ilflj and The staggered structure factor is 
sampled "on the run" as loop updates are constructed 
121 , with implementation of the directed- loop algorithm 
131 to minimize bounce probabilities within a loop. This 
procedure yields a significantly higher sampling efficiency 
than the conventional operator-loop update: close to crit- 
icality (J' > 0.2 J and h < 0.7 J), where the maximum 



number of runs is required, autocorrelation times are re- 
duced by factor of up to 8. Simulations were performed 
for cubic clusters of up to 8 000 spins, and for tempera- 
tures down to r = 0.03J, except at very small J' where 
smaller values are readily obtainable. 

The zero-temperature phase diagram of the coupled 
dimer system is shown in Fig. ^ In the absence of 
a magnetic field the computed quantum critical point, 
separating the spin-disordered, dimer valence bond solid 
(DVBS) phase from the ordered AF phase, occurs at 
(J'/J)c = 0.255(5). Both phases are unpolarized at 
h — 0. Their elementary magnetic excitations have 
the same fundamental difference as for the field-induced 
QPT, except that the coupling- induced ordered phase 
has two gapless Goldstone modes (associated with the 
broken 0(3) symmetry) corresponding to the phase fluc- 
tuations of the staggered moment. The uniform magneti- 
zation varies continuously with applied field from zero in 
the DVBS phase {h < hd) to full polarization at h > hc2- 
In the ordered phase the staggered moment first increases 
with h beyond hd (also for J' > J^) before falling con- 
tinuously to zero at /ic2- 

The values of hc2 obtained from SSE QMC simula- 
tions (triangles in Fig. ^) coincide with the exact re- 
sult, which is classical because quantum fluctuations are 
entirely suppressed in the saturated regime. By con- 
trast, quantum corrections are essential to the determa- 
tion of hci , where a straightforward perturbative expan- 
sion (dashed line) departs from the numerical data (cir- 
cles) for J' > 0.08J. However, the self-consistent bond- 
operator solution (solid line) provides an accurate de- 
scription of the numerical results at all interdimer cou- 
plings, including those close to J^. The quantitative sue- 
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FIG. 3: Critical temperature for field-induced order in 3D 
systems of coupled dimers for a range of values of J'. Sym- 
bols are taken directly from finite-size scaling of SSE QMC 
simulations, except for the T — points which are extracted 
from fits of Tc(h). Solid lines are guides to the eye. 
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FIG. 4; Critical exponents a\ at hci (a) and 02 at hc2 (b), 
shown as a function of the fitting range x (see text). Filled 
circles in (a) are experimental data from Ref. 1^. Dashed 
lines mark the result of the BEC description. 



cess of this theory for KCuCls and TlCuCls may be as- 
cribed to the 3D nature of the ordered phase. 6] 

Turning to finite-temperature properties, the transi- 
tion temperatures for the ordered phase are determined 
by finite-size scaling of the order parameter . Fig. [21 
shows one example of this procedure: in 3D and at the 
critical temperature, at fixed field is expected to obey 
the scaling relation oc L^^~^+'')/^, where the dynam- 
ical exponent z is zero for a thermal transition, [ij] For 
a system with spontaneous breaking of the residual 0(2) 
symmetry the correlation exponent can be taken from 
the classical 3D XY model, for which 77 = 0.0381(2) [l5j |. 
In the vicinity of the critical temperature Tc, the scaled 
m-^ is linear in T |l6l |. and the intersection of the curves 
obtained for different system sizes L determines Tc at all 
fields h G [ft-ci,^c2] for a given interdimer coupling J'. 
Magnification of the intersection region (inset. Fig. \^ 
reveals a set of intercepts within a temperature window 
AT ~ 0.0005J, which determines Tc{h) to very high ac- 
curacy. [T^ 

Fig. O shows the critical temperatures T(.{h) obtained 
for four interdimer coupling ratios, {J' I J) — 0.07, 0.15, 
0.20, and 0.25 which span the range from narrow- 
band, large-gap systems to the wide-band, small-gap 
paradigm. Each coupling ratio corresponds to a vertical 
section in Fig.Q] where the values of hd and hc2 in Fig.|3| 
mark the boundaries of the phase with finite induced 
staggered order; hd vanishes for J' > J'^. The weak 
asymmetry of the Tc{h) curve about the midpoint field, 
{hci+hc2)/2, becomes more pronounced for stronger cou- 
plings due to changes at the extrema of the magnon dis- 
persion relations determining the transitions . The val- 
ues of Tc{h) obtained by the systematic scaling procedure 



confirm, within the error bars, the location of the dips 
in uniform magnetization reported in previous numerical 
studies ^] , and used to construct the experimental phase 
diagram of TlCuClg H. 

We now focus on the scaling properties of Tc{h) close 
to hci and hc2- The scaling exponent a in Eq. (1) is 
obtained by fitting the data shown in Fig. |21 with he and 
a as free parameters. We fit this data within a varying 
window, i.e. for given J'/ J we define 

^ h - hci hc2 - h 
^i(^) = h = 7 T— , (5) 

as normalized measures of the proximity to the lower and 
upper critical fields. Using a minimum of 6 data points 
in a given window xi or X2, the exponents ai and a2 are 
then determined by non-linear curve fitting. The fitted 
exponents are crucially dependent both on the number 
of points and on the exact critical field (Fig. 

The dependence on window size Xi of the exponent 
ai is shown in Fig. 2Ja), along with data extracted from 
experimental reports on TlCuCla A first unambigu- 
ous conclusion from Fig. 0] is that the fitted exponents 
become smaller as the window size is decreased. In par- 
ticular, a large window leads to overestimation of ai, as 
is shown by performing the same varying-window fit to 
the experimental data for Tc{h) [Fig. Ela)]. It is also 
clear that the extrapolated values approach the result 
obtained from the BEC description, Q!bec = 1-5. Higher 
exponents suggested by previous results are seen to be 
a product of the fitting procedure employed, and are in 
fact consistent with the BEC scenario at the transition 
[Fig.^a)]. Most remarkably, the lines for all interdimer 
couplings not only scale to the same limit but do so with 
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the same slope for all x, demonstrating that the critical 
behavior of the system is truly universal for all ratios of 
bandwidth to spin gap. 

Results for a2 are shown in Fig. EJb) , and are again 
consistent with the BEC description as X2 — > 0. Here 
the spread of data points is smaller than for ai{xi), and 
the slope of a2{x2) is less steep. These results may be 
attributed to the less significant role of quantum fluctu- 
ations near hc2, which leads to smaller deviations from a 
mean-field scaling regime. 

It must be emphasized that the universal scaling be- 
havior is obtained only for very narrow windows, meaning 
very low temperatures. Indeed, from Fig. ^ we are un- 
able to determine a temperature below which true scaling 
emerges, and our results indicate that the BEC regime 
is approached only as T ^ 0. At both critical fields, the 
physical origin of departure from universal scaling may be 
found in the approximations inherent to the BEC descrip- 
tion 0- The most important of these are temperature- 
driven renormalization of the quasiparticle effective mass 
{i.e. reduction in curvature at the band minimum) and 
of the effective chemical potential {i.e. band narrowing 
with concomitant increase in spin gap). 

In summary, we have used an unbiased numerical tech- 
nique to study field-induced magnetic ordering transi- 
tions in cubic systems of antiferromagnetically coupled 
dimers. The phase diagram contains two quantum criti- 
cal lines, hci{J' / J) and hc2{J' /J). Both transitions are 
continuous, as is the variation in physical properties of 
the intermediate AF phase, as a function of applied mag- 
netic field. The critical scaling exponents a extracted 
for the transition temperature close to hd and hc2 are 
found to be independent of the interdimer coupling ratio, 
demonstrating that, despite large quantitative differences 
in their niagnon dispersion relations, the physical proper- 
ties of the coupled-dimer systems KCuCla and TlCuCla 
display the same universal scaling behavior. The values 
of a are consistent with a description of the ordering tran- 
sitions as the Bose-Einstein condensation of magnon ex- 
citations, ctBEC = 1-5. However, this regime is obtained 
only in the limit of temperatures very small compared 
to the interdimer coupling, and thus is difhcult to access 
experimentally. Finally, the fits to the exponents were 
shown to depend sensitively on the window size, with 
the consequence that a is easily overestimated. This re- 
sult offers an explanation for the wide range of scaling 
exponents reported in the literature. 
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